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ABSTRACT 

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galax- 
ies at z ~ 10 in a highly overdense region of the universe. These objects could represent the high 
rcdshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs 
at z ~ 6 — 7. To overcome the sampling and resolution problems in cosmological simulations of these 
rare regions, we have used the Constrained Realizations method. Our main result for z ~ 10 shows the 
high-resolution central region of 3.5/i _1 Mpc radius in comoving coordinates being completely domi- 
nated by disk galaxies in the total mass range of i; 10 9 ft. _1 M Q . We have verified that the gaseous and 
stellar disks we identify are robust morphological features, capable of surviving the ongoing merger 
process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to 
negligible numbers. At this redshift, the disks appear to be gas-rich compared to z — 0, and the dark 
matter halos baryon-rich, by a factor of ~ 2 — 3 above the average fraction of baryons in the universe. 
The dominance of disk galaxies in the high density peaks during the epoch of reionization is contrary 
to the morphology-density trend observed at low redshifts. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: halos — galaxies: interactions 
— galaxies: kinematics and dynamics — galaxies: star formation 



1. INTRODUCTION 

Deep sky surveys are providing large sample of color- 
selected galaxies out to z ~ 6 (e.g., Bouwens et al. 2007) 
and a growing number of spectroscopically confirmed ob- 
jects (e.g., Fan et al. 2004; Malhotra et al. 2005; Stark et 
al. 2011). With the availability of the Wide Field Cam- 
era 3 (WFC3) on the Hubble Space Telescope (HST) 
more than 100 galaxies have been detected at z ~ 6 — 8.5 
(e.g., Bouwens et al. 2010; McLure et al. 2010; Oesch 
et al. 2010). These high-z galaxies are expected to be 
significantly clustered in the high density environments 
(e.g., Overzier et al. 2006). 

Recent studies of high-redshift QSOs have indicated 
that they reside in rare and highly overdense regions at 
z ~ 6, with the comoving space density of ~ (2.2 ± 
0.73) h 3 Gpc" 3 (e.g., Fan et al. 2004). These dense en- 
vironments, inhabiting ~ 1 Gyr old Universe, can be 
precursors of massive clusters of galaxies at z = (e.g., 
Fan et al. 2004; Springel et al. 2005; Li et al. 2007), but 
not necessarily the most massive clusters (e.g., Trenti et 
al. 2008; Romano-Diaz et al. 2011). Chemical compo- 
sition inferred from QSO spectra and additional correla- 
tions with the cosmic star formation history indicate a 
possible link between galaxy evolution and the QSOs — 
being metal-rich they plausibly are located at the cen- 
ters of massive galaxies (Graham et al. 2001; Barth et 
al. 2003; Bunker et al. 2004). 

Galaxy evolution in these overdensities can differ from 
that in the field, both quantitatively and qualitatively, 
due to earlier formation times, an absence or a decreased 
intensity of the background UV field, and a possible ra- 
diative and mechanical feedback from the QSOs. The 
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reionization of the Universe could start with these re- 
gions as early as z ~ 15 and extend to z ~ 6 (e.g., Becker 
et al. 2001; Barkana 2002; Cen & McDonald 2002). But 
the population of galaxies responsible for re-ionization 
of the Universe is not yet found. It is also not clear at 
present what fraction of overdense regions and to what 
extent is affected by QSOs and background UV. 

The most overdense peaks are hard to sample with 
cosmological simulations at high resolution. Recent de- 
tections of candidate galaxies and proto-clusters at z ~ 
6 — 10, have inspired us to study properties of galaxies 
at these redshifts (e.g., Trenti et al. 2011a, b; Bouwens 
et al. 2011; Oesch et al. 2011; Capak et al. 2011), and 
this is first in the series of papers on this issue. Cir- 
cumventing the difficulty of sampling, we resort to the 
Constrained Realizations (CR) method (e.g., Hoffman & 
Ribak 1991), a fully self-consistent method which allows 
to simulate such rare regions at highest resolution, with- 
out loss of generality. 

Current efforts in modeling galaxy evolution aim at re- 
producing the observed morphology at z = (e.g., Gov- 
ernato et al. 2004; Robertson et al. 2004; Tasker & 
Bryan 2008), but also analyze the relevant physical pro- 
cesses at high z (e.g., Gnedin et al. 2009; Romano-Diaz 
et al. 2009; Gnedin & Kravtsov 2010; Greif et al. 2010; 
Wise et al. 2010; Pawlik et al. 2011). 

To probe the dark matter (DM) backbone of galaxy 
formation at z ~ 6 — 10 in highly overdense regions, 
Romano-Diaz et al. (2011) have performed a set of care- 
fully designed DM simulations with CRs. Constructed 
halo Mass Functions have shown that on average, the 
QSO hosts have enhanced the halo formation in their 
vicinity due to pure gravitational effects, when compared 
with an average region in the Universe. Here we extend 
our study to include treatment of the baryonic compo- 
nent, and analyze the properties of galaxies at z ~ 10. 
In this Letter, we focus on the most general properties 
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of galaxy population, such as morphology and baryon 
fraction. In the long run, our goal is to quantify the en- 
vironment of high-z QSOs, understanding the similarities 
and peculiarities of galaxy evolution in the overdensest 
regions at z ~ 6 — 10. While current simulation does not 
include the QSO feedback, and for z ~ 10 neglects the 
UV background, it would be interesting to account for 
these effects in the future. 

This Letter is structured as follows: section 2 summa- 
rizes the technical details of our simulation setup and 
its initial conditions; section 3 analyzes obtained galaxy 
morphology at z ~ 10, and the last section discusses the 
obtained results. 

2. NUMERICS AND INITIAL CONDITIONS 

We use the modified version of the tree-particle- 
mesh Smoothed Particle Hydrodynamics (SPH) code 
GADGET-3 originally described in Springel (2005) , in its 
conservative entropy formulation (Springel & Hernquist 
2002). Our conventional code includes radiative cool- 
ing by H, He, and metals (Choi & Nagamine 2009), star 
formation, its feedback, a phenomenological model for 
galactic winds, and sub-resolution model of multiphase 
interstellar medium, ISM (Springel & Hernquist 2003). 
Since we only focus on high redshift before the full reion- 
ization and do not implement on-the-fly radiative trans- 
fer of ionizing photons, the UV background is not in- 
cluded in our simulation. In the multiphase ISM model, 
star forming SPH particles contains the cold phase that 
forms stars and the hot phase that results from a super- 
nova (SN) heating. The cold phase contributes to the gas 
mass, and the hot phase contributes to gas pressure. For 
the star formation, we use the "Pressure model" which 
reduces the high-z star formation rate (Choi & Nagamine 
2010) relative to the previous model by Springel & Hern- 
quist (2003). The star formation is triggered when the 
gas density is above the threshold tih.sf = 0.6 cm' 3 . 

The initial conditions have been generated using the 
CR method (Hoffman & Ribak 1991; see also Romano- 
Diaz et al. 2007) and are similar to those used by 
Romano-Diaz et al. (2011, see the Appendix there). 
The constraints were imposed onto a grid of 1,024 3 
within cubic box of 20h~ 1 Mpc to create a DM halo seed 
of 10 12 h~ 1 MQ collapsing by z ~ 6, according to the 
top-hat model. We assume the ACDM cosmology with 
WMAP5 parameters (Dunkley et al. 2009), fi m = 0.28, 
n A = 0.72, ft h = 0.045, and ft = 0.701, where ft is the 
Hubble constant in units of 100 km s _1 Mpc -1 . The vari- 
ance us = 0.817 of the density field convolved with the 
top hat window of radius 8ft.' 1 Mpc -1 was used to nor- 
malize the power spectrum. 

The evolution has been followed from z = 199. The 
growth of structure has been substantially accelerated 
by resulting high overdensity, ~ 20, within the central 
region. The simulation runs in comoving coordinates and 
with vacuum boundary conditions. We have applied the 
zoom-in technique with 3 levels of refinement to increase 
the numerical resolution. The inner region has a radius 
3.5ft' 1 Mpc, having an effective resolution of 1,024 3 in 
DM and SPH particles. Within this region, we obtain 
the particle mass of 4.66 x 10 5 M Q (DM), 1.11 x 10 5 M Q 
(gas) and 5.55 x 10 4 M Q (stars). Gravitational softening 



stein & Hut 1998) to identify DM halos, and have defined 
halo virial masses and radii in the context of the spherical 
top-hat collapse model, M v ; r = 4/37rA(z)/9(z)-R 3 ir , where 
p is the cosmic background density and A(z) (Bryan & 
Norman 1998) is the critical overdensity at virialization. 
The halo and disk shapes have been found using the 
method described in Romano-Diaz et al. (2009) follow- 
ing Heller et al. (2007). 
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We have used the halo finding algorithm HOP (Eisen- 



FlG. 1. — Example of the phase diagram T e g — n, effective gas 
temperature vs density, for the gas within the most massive DM 
halo at z = 10.2. Density is normalized by jih.sf = 0.6 cm -3 
— the critical density for the onset of star formation (see text). 
The colors are: disk gas particles (red), halo gas particles (blue), 
and wind particles (black). The blue-red boundary represents the 
empirically determined transition between the disk and halo gas 
(see text). This snapshot corresponds to a post-merger evolution 
with a spike in star formation rate. 



3. DISK FRACTION 

By z ~ 10, nearly 500 halos above lO 8 ft _1 M0 have 
been detected and cataloged. To analyze the baryons 
within the DM halos, they have been divided into three 
groups: (1) gas; (2) stars; and (3) gas which currently ex- 
periences feedback from stellar evolution and is found on 
ballistic trajectories — this is the wind. To identify and 
quantify gas properties, we have constructed the T c g — n 
phase diagram (e.g., Fig. 1). The gas tends to cluster in 
two main regions: the top ones (blue and red, see below) 
are the disk/halo particles which can experience the star 
formation and whose effective temperature, T e g, is the 
result of a shock heating by the SN. The second region 
(black) consists of particles that comprise the wind. All 
gas can be cooled/heated via expansion/compression. 

We do not know apriori whether a halo contains a gas 
disk, and, therefore, identify its presence, based on axial 
ratios determined from the moments of inertia (Heller et 
al. 2007; Romano-Diaz et al. 2009). Major and minor 
equatorial axes are defined as a and 6, and c (disk semi- 
thickness) along rotation axis of the disk. Specifically, we 
focus on the ratio of the minor-to-major axis, c/a — its 
actual values are discussed below. If, based on the c/a 
ratio, we decide that a particular halo does not contain 
a gas disk, we call the gaseous entity irregular. We have 
empirically determined that disk and halo gas can be 
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Fig. 2. — Example of a deprojected face- and edge-on strongly barred disk embedded in the most massive halo at 2 ~ 10.2, stars (top) 
and gas (bottom). The frames are 5x5 h~ 1 kpc. Left frames: bar along the line-of-sight; right frames: bar normal to the line-of-sight. The 
total mass is 

JWtot ~ 1-11 x lO 10 ^- 1 M Q , the total disk mass is ~ 2.9 X 10 9 /i 1 Mq, The gas disk shown corresponds to the red colored 
zone in Fig. 1. 



separated by a density cut at n/nn.SF ~ 100- The red 
color in Fig. 1 corresponds to the disk gas and the blue 
color to the halo gas. We have verified this separation 
randomly for a large sample of DM halos. Of course, in 
the absence of a disk this separation is not important, 
but we still keep it for the sake of simplicity. The disk 
stars have been separated from halo stars based on the 
alignment of their angular momenta (e.g., Romano-Diaz 
et al. 2009). 

After isolating the disk gas, we have deprojected its 
spatial distribution in each halo onto three fundamental 
planes: face-on and edge-on (see an example in Fig. 2 
for the most massive halo) . The disks are clearly visible, 
have sizes i?d ^ 3/i _1 kpc and are embedded in halos 
of i? v ; r <; 7/i _1 kpc. The face-on disks show large-scale 
bars, sometimes comparable to the disk size, and their 
gas layer full thickness is ~ 1 — 2/i _1 kpc, and is pressure- 
supported. 

Figure 3a shows the disk gas and stars, and DM frac- 
tions as functions of the total halo masses. The cutoff 
at logM to t ~ 8.0 approximately corresponds to ~ 350 
particles within i? v ; r and we consider this to be reliable 
in determining halo (e.g., Trenti et al. 2010) and disk 
shapes. Two types of behavior can be observed in Fig. 3a. 
The DM fraction drops and gas fraction increases to 
logMtot ~ 8.5, and the stellar fraction increases up to 
logMtot ~ 9.5. The DM and stellar fractions flatten 
thereafter, while the gas fraction declines. 

The stellar disks are generally less massive than the 
gaseous ones for most of the mass range, logM to t Si 9.5 
(Figs. 2b, 3a). The overall trend is a strong increase 
in the gas fraction in disks with the total mass, up to 



logAftot ~ 8.5, to ~ 1/3, and a decrease thereafter. The 
steep decline in the gas fraction below logMtot ~ 8.5 is 
related to the inability of the halo to retain all of the gas 
because of the SN feedback and low virial temperatures. 
Below logMtot ~ 8.0 (not shown here), the gas fraction 
in the halos is negligible because of the same effect. For 
the most massive halos, there appears to be roughly equal 
fraction of disk gas and stars, ~ 50%. 

The shapes of the gaseous disks show interesting trends 
(Fig. 3b). Overall disks appear mildly non-axisymmetric, 
0.8 < b/a <S 1, with a slight increase in the dispersion 
of b/a toward smaller M to t- On the other hand, there 
appears to be a strong trend toward increased flatness, 
c/a, with increasing mass, down to c/a ~ 0.1. Below 
logMtot ~ 8.3, the gaseous disks appear prolate and 
their rotation axis can be distinguished only based on 
the angular momentum. We also note that parent halos 
of galaxies discussed here are quite triaxial with the range 
of b/a ~ 0.6 — 1 and c/a ~ 0.2 — 0.8 and therefore can 
influence the shape of the disk they harbor (Berentzen 
& Shlosman 2006). 

The stellar disks (e.g., Fig. 2b) exhibit a high den- 
sity midplane component of newly formed stars, with 
the mass of ~ 8 x 10 8 /i _1 M© for the most massive halo 
at z = 10.2, and a lower density spheroidal stellar com- 
ponent, of ~ 1.2 x 10 9 /i _1 Mq, when integrated over the 
extended halo. Note that the star formation we measure 
is not limited to the disk regions only, but engulfs the en- 
tire halos because the gas density there appears above 
the threshold density nn.SF- So the spheroidal compo- 
nent in our galaxies has its origin in the halo and was not 
formed as a result of the secular or dynamical thickening 
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Fig. 3. — Properties of halos and galaxies at z ~ 10.2 within 
the high-resolution region of 3.5/i -1 Mpc in radius, with binning 
of 0.25 in log (Aftot Mq) ^ 9.0 and 0.1 dex for smaller masses, in 
2 X 1,024 3 simulation, (a) disk gas, disk stars and DM fractions 
within halos as a function of the halo total mass (species: DM — 
black, stars — red, gas — blue). The arrow shows the WMAP 
baryon fraction, f2 D /f2 m ; (b) gas disks axial ratios b/a (blue) and 
c/a (red); (c) disk fraction, /j, as a function of the halo total mass, 
Aftot- Shown are two defining criteria for the gas disk, c/a £ 0.5 
(blue) and c/a 0.7 (red). Halos in the process of merging at this 
z have been omitted, overall 4 objects within this mass range — 
all of them had more than one disk per halo. The control 2 X 512 3 
simulation gave nearly identical results for Fig. 3c. 



of the stellar disk, or from the heating by minor mergers. 
It is not clear what is the fate of this spheroidal stellar 
component. Does it represents low metallicity halo stars 
at low redshifts? Following the evolution of the relative 
masses of disk and halo stars, we only comment that the 
spheroids have formed early, by z ~ 12, and we observe 
a steadily increasing disk masses, compared to the halo 
stars, with time. 

We are now in a position to estimate the disk frac- 
tion, /<j, within the high resolution box of the simula- 
tion. To show the robustness of our result, we adopt 
two definitions of the gaseous disks, these of c/a 0.5 
and c/a SS 0.7. The first and more conservative def- 
inition exhibits a sharp decline below logM to t ~ 9, 
while the second one shows a very similar decline below 
logM tot ~ 8.6. Hence, both agree that the disk fraction 
is 100% for the more massive halos (Fig. 3c). There is 
some indication that the disk fraction levels off again for 



very small halos, at the level of /<j ~ 0.05 — 0.2. 

We have checked the sensitivity of /d to our definition 
of the disk gas threshold (§3). Varying the threshold by 
±10% does not change /d for logMtot ^ 9, and changes 
it negligibly for smaller masses. Furthermore, we have 
tested /d in lower resolution 2 x 512 3 simulation and 
find it nearly identical — a sign of numerical convergence. 
Finally, we note that the blue and red point scatter is 
uniform across the mass range in Fig. 3b — a sign that 
the resolution effects do not compromise the results. In 
summary, the disk fraction is /d ~ 1 for for massive 
halos in our simulation, and exhibits a sharp decline by 
an order of magnitude below logM to t ~ 8.5 — 9.0. 

The circular velocities for massive halos are con- 
tributed by baryons and DM nearly equally. Only for 
the most massive halo, we see maximal v c J> 100 km s -1 . 
For the 15 most massive halos, we have instead maxi- 
mal v c 50kms _1 . Overall, the total mass in a halo 
roughly correlates with the maximal v c , meaning that 
the halo concentration is a monotonic function of M to t. 

4. DISCUSSION 

Using a high-resolution cosmological numerical simula- 
tion, we have analyzed the evolution of galaxies at z ~ 10 
in a highly overdense region of the universe. These ob- 
jects could represent high redshift galaxies recently ob- 
served by the WFC3/HST, and be possible precursors of 
QSOs at z ~ 6 — 7. For this purpose we have employed 
high-resolution simulations of 3.5/i _1 Mpc region in ra- 
dius involving a multiphase ISM and feedback from stel- 
lar evolution. Halos with M to t ~ 10 8 /i _1 M are resolved 
with ~ 350 particles. We find that the fraction of disk 
galaxies is unity for total masses above ~ 10 9 /i _1 Mq, 
with a sharp decline in /d for less massive halos. We 
also find that the baryon fraction (with respect to total 
mass) peaks around logMtot ~ 8.5, slowly declines to- 
wards higher masses and sharply declines toward smaller 
masses. So overall, above the resolved total mass of 
10 8 /i _1 M , halos are baryon-rich at these redshifts with 
respect to the WMAP f2 D /^m ratio, in a broad agree- 
ment with Harford et al. (2008). More massive halos 
appear with about equal gas and stellar fractions, while 
less massive ones have higher gas fractions. 

The prevalence of disk galaxies at high redshifts is not a 
trivial result because of ongoing merging between galax- 
ies which contributes to strong perturbations in disk mor- 
phology, and has been tracked by us. Four objects in the 
process of merging have been omitted from our statistics 
for /d- All the mergers we observe are between disks. De- 
spite this process, we find that gas-rich halos are capable 
of reforming the disks in a short time. This effect has 
been demonstrated previously by Springel & Hernquist 
(2005) and we confirm it here — the mergers we have ob- 
served do not destroy the disks. Taken at the face value, 
the dominant population of disk galaxies at overdense 
peaks at z ~ 10 appears to contradict the morphology- 
density relation observed at low z (e.g., Dressier 1980). 
The most straightforward explanation for this is that var- 
ious processes leading to deficiency of disks in the density 
peaks at z = (e.g., dry mergers, gas ablation) have had 
not enough time to operate at z ~ 10. A number of 
processes at high z contribute to galaxies being gas-rich, 
among them the ongoing gas accretion which joins the 
disks in their equatorial planes. In tandem with promi- 
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nent (gas-rich) bars, these disks can serve as a fuel reser- 
voir for growing supermassive black holes (SBHs) already 
at this epoch. 

While it has been claimed that the feedback from such 
SBHs (neglected here) can terminate the accretion pro- 
cess and drive the gas out of the galaxy (e.g., Di Matteo 
et al. 2005), it is by no means a solved problem. It is 
also not clear when the morphology-density correlation 
observed here is reversed to its low-redshift counterpart. 

Our simulation has demonstrated that well resolved 
galaxies hosted by logAf to t ^ 9 mass halos possess well 
defined gaseous and stellar disks. Moreover, by observing 
the evolution of these galaxies through major mergers, we 
have verified that the disks are also robust morphologi- 
cal features and survive this process. We do agree with 
Greif et al. (2010) and Wise et al. (2010) simulations of 
smaller halos that the stellar feedback becomes increas- 
ingly dominant when smaller, ~ 10 s hr 1 Mq halos are an- 
alyzed, and with Pawlik et al. (2011) who have obtained 
the disk formation in higher mass halos. We interpret 
the increased irregular morphology below logM to t ~ 9.0 
as the effect of a stellar (SN) feedback on the disk gas. 
At this mass range, the morphology becomes chaotic and 
the gas fraction declines sharply. The observed decreas- 
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ing disk fraction in low mass halos appears to be real. 
However, the analysis of the feedback process to quan- 
tify this phenomenon is outside the scope of this Letter. 

The dominant population of the disk galaxies at z ~ 10 
in the overdense regions simulated in this work appear 
consistent with the presence of extended morphologies 
in bright z ~ 7 galaxies observed by HST (Oesch et al. 
2010). However, a direct comparison of our predictions 
will be possible only with next generation telescopes such 
as JWST or 30m-class observatories from the ground, 
which will have sufficient sensitivity and angular resolu- 
tion to investigate the morphology of z <; 8 galaxies. 
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